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ABSTRACT 

There is a significant nonlinear correlation between the Eddington ratio (iboi/^Edd) 
and the Eddington-scaled kinetic power (ikin/^Edd) of jets in low luminosity active galac- 
tic nuclei (AGNs) (Merloni & Heinz). It is believed that these low luminosity AGNs contain 
advection dominated accretion flows (ADAFs). We adopt the ADAF model developed by Li 
& Cao, in which the global dynamics of ADAFs with magnetically driven outflows is de- 
rived numerically, to explore the relation between bolometric luminosity and kinetic power of 
jets. We find that the observed relation, Lkin/^Edd oc (-Lboi/iEdd)""*^, can be well repro- 
duced by the model calculations with reasonable parameters for ADAFs with magnetically 
driven outflows. Our model calculations is always consistent with the slope of the correlation 
independent of the values of the parameters adopted. Compared with the observations, our 
results show that over 60% of the accreted gas at the outer radius escapes from the accretion 
disc in a wind before the gas falls into the black holes. The observed correlation between 
Eddington-scaled kinetic power and Bondi power can also be qualitatively reproduced by our 
model calculations. Our results show that the mechanical efficiency e (Lkin = e^^bondic^) 
varies from 10^^ 
simulations. 



10 , which is roughly consistent with that required in AGN feedback 
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1 INTRODUCTION 

It was found that the mass of the black holes in the centre of galax- 
ies is t ightly correlate d with galaxy properties (e.g.. |Tremaine et ahl 
I2OO2I : iBandara et al.j .2009.) . The co-evolution of massive black 
holes and galaxies is usu ally ascribed to the so-called AGN feed- 
back jSilk&Reeslll998l) . There are two ways of the AGN feed- 
back on the ambient gas, i.e., radiative feedback and mechan- 
ical feedback, which correspond to the case s of high and low 
Eddington-scaled accretion rate respectively jpabian et alj|200ll ; 
ICrotonetalJl2006l) . The mechanical feedback is in the form of 
powerful outflows/jets produced in the most central small region, 
which can spread to a large region. The X-ray cavities and bub- 
bles in the galaxies and clusters are believed to be blown up 
through the interaction between the outflows/jets and the intra- 
cluster medium (ICM), which are carefully studied owing to the 
high resol ution observations of Chandra and XMM-Newton (e.g., 
Birzan et a l. 2004; Allen et a l. 2004, 2006; Rafferty et al. 20q3; 



containing the data of mechanical power, the black hole mass and 
the nuclear luminosity both in the radio (5 GHZ) band and in the 
2-10 kev band. They found a strong correlation between Eddington- 
scaled kinetic power and bolometric luminosity. A similar correla- 
tion between Eddington-scaled kinetic power and Bondi power was 
also present for this sample (see also[Allen et al. 2006). 

ADAFs are optically thin, geometrically thick and very hot 
as most releasing gravitational energy being stored in the gas in- 
stead of radiating away, which are believed to be present in most 
of low-luminosity AGNs (e.g.. lNaravan & Yilll994ll995h . Due to 
their high temperature and positive Bemouli parameter, ADAFs ar e 
prone to be associated by outflows jBlandford & Begelman|[T999h , 
which are consistent with numerical simulations and observations 



(e.g..lStone & Pringlel200ll ; lGammie et aljl999l ; l\\ian et al.l2003h . 



Merloni & Heinzli2007l). S earch t he literature for the data of X-ray 
cavities. iMerloni & Heinj ( l2007h complied a sample of 15 objects 
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iMerloni & Heinj " 2007 ) adopted a simple sel f- similar model for 
ADAFs coupled with outflows developed by iKuncic & Bicknell 
( I2OO4I) to explore the correlations between Eddington-scaled ki- 
netic power, bolometric luminosity and Bondi power. They found 
that the correlati ons can be qualitatively reproduced by the model 
calculations (see lMerloni & Heinzll2007l for the details). 

The outflows may probably be accelerated from the accre- 
tion flow by the larg e-scale ordered magnetic fields threading the 
accretion flow (e.g., iBlandford & Pavnell 19821 ; Icaoll2002l ; ISpruitl 
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I2OO8I : iLi, Gan.&Wan3l2010l . iLi & Cao|[2009l) . iLi & Cad | |2009|) 
studied the global dynamics of ADAFs with magnetically driven 
outflows. In their calculations, the strength of the large-scale mag- 
netic fields is assumed to scale with the gas pressure of the ac- 
cretion flow. The global structure of the accretion flow together 
with the mass loss rate in the outflows are derived simultaneously, 
and then the kinetic power of the outflow is available. In this pa- 

i ier, we a dopt the model of ADAFs with magnetically outflows 
Li & Caoii2009 ) to explore the above-mentioned correlations given 
bv lMerloni & Heinj l2007h . 



the ADAF can be described by one-dimensional hydrodynamical 
equations jNaravan & YillT99a ). The torque exerted on the accre- 



2 MODEL 

In this work, we adopt the model of ADAFs with magnetically 
driven outflo ws/jets surroundin g a non-rotating black hole, which is 
suggested by |Li & Cad l l2009t) . In this model, the physical quanti- 
ties of the ADAF are integrated in vertical direction, and therefore 
F can be described by 
jNaravan &Yill 19951) . 
tion flow due to the magnetically driven outflow. In principle, the 
connection between the accretion disc and the outflow can be ex- 
plored if the magnetic fields threading the accretion disc is known 
and the vertical structure of the d isc is provided as the bound- 
ary condition of the outflow (e.g., ICao&SpruiJll994l) . For sim- 
plicity, we assume that the pressure of the magnetic fields at the 
disc surface is proportional to the total pressure of the accretion 
disc in this work. The mass loss rate in the outflow is evaluated 
by assuming that the strength of the magnetic fields threading the 
disc to be roughly self-si milar above the disc as that adopted in 
iBlandford & Pavnd ( fl982l) . and the effects of the outflow are then 
properly included in the one-dimensional hydrodynamical equa- 
tions for the ADAF. Although this treatment is rather simplified, 
we believe that it should be sufficient good for modeling the ob- 
served statistic correlations in radio galaxies. 

We summarize this model briefly in this section. The basic 
equations, i.e., the continuity equation, radial momentum equation, 
angular momentum equation and the energy equations are given as 
follows. 

In the presence of magnetically driven outflows, we need to 
add some additional terms related with ou tflows to the dynamical 
equat ions for a classic ADAF (see ICao & Spruit 2002; Li & Caol 
I2OO9I for the details). The continuity equation describing such an 
ADAF+outflow system is 

-4(27ri?EwR) + 47ri?mw = 0, (1) 

where all the physical quantities denote their common meanings, 
rhw is the mass loss rate from unit surface area of accretion flow 
due to the outflow. 

The radial momentum equation is 

"R-jB-^-R^ - ^^k) + - 5m = 0, (2) 
dti p au 

where Q is the angular velocity of the accretion flow at R, the radial 
magnetic force Qm = -BrBz/SvrE (B^ and B^, are the radial and 
vertical components of the magnetic fields at the disc surface). 
The angular momentum equation reads 



VR- 



dR 



(3) 



where T,„ is the magnetic torque exerted on the unit sur- 
face area of the disc due to the presence of the outflows 



jLubow, Papaloizou. & Pringldll994bl) . A fraction of the energy 
and angular momentum of the accretion flow is carried away by the 
outflows, which is equivalent to a torque exerted on the accretion 
flow. We can calculate the magnetic torque Tm with energy con- 
servation law, provided the power of the outflow from unit surface 
area of the accretion disc is known, which leads to 



1 



(4) 



were /kin and rhw are the kinetic power and the mass loss rate of 
the outflow from unit surface area of the accretion disc, Sc and ei 
are the specific internal energy of electrons and ions respectively. 
The main difference between our cal culation of Tin and that in 
iLubow. Papaloizou. & Pringld l ll994bl) 's work is that the internal 
energy of the gas is properly included in our calculation, while 
theirs is only for cold outflows. 

The energy equations for ions and electrons are given by 
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0, 



(5) 



2mw£i 
2H 



0, 



(6) 



and 

,dei Pi dp ^ ^,1 

P^'^^dR-^PdM^-^'^'^'^ + 

where, the parameter 5 describes the fraction of the viscously dissi- 
pated energy that goes directly into electrons in the accretion flow, 
= —aPRdQ,/dR is the energy dissipation rate per unit vol- 
ume, and the radiative cooling rate con sists of synchrot ron, 
bremsstrahlung, and Compton coolings (see, iManmotollioOol for 
details). The last terms in Eqs. l[5) and ^ represent the internal en- 
ergy of the gas in the accretion flow carried away in the outflows. 
The effects of bulk kinetic part of the outflows are included in the 
angular momentum equation l[3]l, and the structure of the ADAF 
is altered, which leads to changes of in Eqs. lO and The 
energy transfer rate gio from ions to electrons through Coulomb 
collisions is 



3 rric 
2 m, 



-TlcTliCrTC 



K2{l/ec)K2{l/ei 



■In A 



Ai 



+ 2Ko 



(7) 



where the Coulomb logarithm InA — 20 dStepnev & Guilber3 
Il983h . 

For simplicity, the large-scale magnetic field lines are assumed 
to thread the accretion disc, and the strength of the magnetic fields 
far fro m the disc surface along the field line to be roughly self- 
similar jLubow. Papaloizou. & Pringldll994bh , 



Bp{R) - Bpd{Rd 



(- 

\Rd 



(8) 



where Bp^ is the strength of the poloidal component of the field at 
the disc surface, Bp{R) is the field strength at R along the field line 
threading the disc surface at _Rd- The self-similar index ^ 1 is re- 
quired, because the configuration of any magnetic fields of the disc 
being able to accelerate outflows should have an expanding shape 
above the disc. In our model calculations, we have not adopted 
a detailed configuration of magnetic fields threading the disc, in- 
stead, we use the parameter to describe how the poloidal field 
str ength decline along the field line for si mplicity as that adopted 
by Lubow. Papaloizou. & Pringle 1 1994bl) . In the calculations of 
iLubow. Papaloizou. & Pringld ( 1 1994bir ^ = 4 is adopted . 

For a relativistic jet accelerated by the m agnetic field of th e 
disc, the Alfven velocity at the Alfven point is jCamenzin 3ll986l) 
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VA = 



(47rpA7f 



(9) 



where Bp and pA are the poloidal field strength and the density of 
the outflow/jet at Alfven point, and is the Lorentz factor of the 
bulk motion of the outflows/jets at the Alfven point. 

The mass and magnetic flux conservation along the field line 
requires 



?Tt„ ^ PAVA 

Bpd 



-Dp 



(10) 



The outflow/jet can be magnetically accelerated over the 
Alfven p oint all along until the modi fie d fast magne t osonic 
surface iU. Chiueh. & Begelmad 1 19921 ; ICao & Spruil 1 1994 
IVlahakis & Konigll I200I l2004h . and the bulk velocity of the 
outflow/jet at the Alfven point has reached a significant frac- 
tion of its terminal value (Spruit 2008). In this work, we have 
not derived an outflow solution passing smoothly through the 
Alfven and slow/fas t magnetosonic points as those done in the 
previ ous works (e.g..lLi. Chiueh & Begelmaniri992l :lc ao & Spruill 



1 1994 IVlahakis & Konigll I200I l2004h ." The Lorentz factor of the 
outflow/jet is 



1_ Ifl 



(11) 



Combining equations dst- dl It . the mass loss rate in the out- 
flow/jet from the unit surface area of the disc is 



47rc 



7j 



(^a2 _ ;^^(l+c)/2 



(12) 



The origin of the ordered magnetic fields threading the 
disc is still unclear. It was suggested that the magnetic fields 
can be generated through dynamo processes in the disc (e 



IShakura & Sunvaevll 19731 ; I Armitagdl 19981 ; iRomanova etai][T99: 
or the large-scale external magnetic fields are transported inward 
by the accretion flow (e.g..lBisnovatvi-Ko gan & Ruzmaikin" 19761; 
iLubow. Papaloizou. & Pringld Il994al ; ISpruit & Uzdenskv 2005). 
For simplicity, we assume the magnetic pressure is proportional to 
the gas pressure in the accretion flow: 



P 



1~I3 



(13) 



where /3 is the ratio of the gas pressure to the total pressure, B is 
the strength of the magnetic fields in the accretion flow. 

The kinetic power of the outflow/jet is extracted from the 
ADAF with magnetic fields. The kinetic power of the outflow at 
Alfven point can be calculated when its mass loss rate is derived, 



Rout 



/kin47r7?d7?, 



(14) 



where Rin and -Rout are the inner radius and outer radius of the 
accretion disc, /kin is the kinetic power of the outflow accelerated 
from the unit surface area of the accretion disc. The cooling of the 
gases during the acceleration is neglected in Eq. ( I14t . The outflow 
can still be accelerated beyond the Alfven point in the outflow (e.g., 
iLi, Chiueh, & Begelmanlfl992h . The kinetic power of the outflow 
from unit surface area of the accretion flow is: 



/kin = (7j - l)mwC^ + 7j niw(ei + £c) + S, 



(15) 



assumed finally to convert to kinetic power of the outflow. For sim- 
plicity, we estimate the Poynting flux S at the Alfven point roughly 
with S ~ -B^^UA/(47r) ( ISpruill2008h . We believe that it should be 
sufficient good for modeling the observed statistic correlations in 
radio galaxies. 

The bolometric luminosity of the ADAF can be calculated as 

n-Rout 

Lbol = / q^AnRHAR, (16) 



where q is the radiative cooling rate cons isting of synchrotron , 
bremsstrahlung , and Compton coolings jNaravan & Yil 1 19951 ; 
iManmotdlioog) . 



3 RESULTS 

Integrating these equations (|Tll-(|6]l from the outer boundary of the 
flow at Tiout inwards the black hole, we can obtain the global struc- 
ture of the accretion flow passing the sonic point smoothly to the 
black hole horizon. In all our calculations, we adopt the black hole 
mass M = 8.8 x lO ^Mp), which is the mean black hole mass of 
the samples given in iMerloni & Heinj j2007h . The conventional 
values of disc parameters are adopte d as: a = 0.1 , 5 = 0.1, 
= 5000-Rg and «a = "uk (see IlI & Caoll2009l for the de- 
tailed discussion). 

As described in Sect.[2l the structure of the accretion flow and 
the outflows is calculated, and then the kinetic power of outflows 
and the bolometric luminosity of the accretion flow are available. 
Changing the mass accretion rate, the relation between Eddington- 
scaled bolometric luminosity A (A = -i/boi/^Edd) and kinetic 
power of the outflows is plotted in Fig. [T] for different values of 
C with P = 0.92, i.e., = 0.087Pgas, in which the results are 
compare d with the correlation betw een Lboi/^Edd and LkinZ-f/Edd 
given by iMerloni & Heinj ^2mh . The Eddington-scaled kinetic 
power as functions of Eddington-scaled bolometric luminosity of 
the accretion disc is plotted in Fig.|2]for different values of /3 (i^ = 2 
is adopted). The observed correlation can be well reproduced by 
our model calculations by tuning the values of two parameters, /? 
and In Fig. [3] the best fitted model parameters in C,-I3 space are 
plotted. We plot the mass accretion rate as functions of radius for 
different values of in Fig.|4] 

The Bondi power is defined as l lMerloni & Heinzll2007h : 



-Pbondi — O.lMbondiC 



(17) 



where A/bondi is the Bondi accret ion rate. The Bondi acc retion 
rates of the sources in the sample of IMerloni & Heiiiz |2007l) were 
estimated with the X-ray observations jAllen et al.ll2006l) . In this 
work, we simply adopt Mq — Mbondi, where Mo is the accretion 
rate at the outer radius. We can therefore plot the dependence of the 
kinetic powers as functions of Bondi power in Fig. |5] for different 
values of /? and ^. 



4 DISCUSSION 



where the Poynting flux S at the Alfven point in the outflow is 



IMerloni & Heinz! l l2007l) 's sample are limited to low luminosity ra- 
dio galaxies, which are supposed to be accreting at low rates, prob- 
ably through a radiatively inefficient ADAF. In t his wo rk, we show 
that the ADAF model developed by iLi & Cad ( l2009h can explain 
the correlation between Ed dington-scaled kinetic po wer and bolo- 
metric luminosity found bv lMerloni & Heina ( l2007h . 

In Fig.[T] we show the kinetic power as functions of bolometric 
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Figure 1. The Eddington-scaled kinetic power as functions of Eddington- 
scaled bolometric luminosity (A = ^bol/^Edd) for different val- 
ues of (. The bo ld solid hne is the b est fitted result of on the 
correlation givn in iMerloni & Hein j i2007h . i.e., log(L]jin/LEdd) = 
0.491og(I/bol/-^'Edd) ~ 0.78. The dashed line, solid line and dotted line 
are for f = 1, f = 2 and C = 4 respectively. 
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Figure 2. The Eddington-scaled kinetic power as functions of Eddington- 
scaled bolometric luminosity (A = ibol/^Edd) for different values of 
/3. The bold sohd line i s the best fitted result of on the correlation givn in 
iMerloni & Heinj j2007h . The dashed line, solid line and dotted line are for 
/3 = 0.9, /3 = 0.92 and /3 = 0.94 respectively. 
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Figure 3. The best fitted results to the observed correlation in the space. 
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Figure 4. The mass accretion rate as functions of radius for different values 
of C^, where /3 = 0.92 and the accretion rate at outer radius mo = 10^''' 
(mo = Ma/MEdd) are adopted. 
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Figure 5. The Eddington-scaled kinetic power as a function of Eddington- 
scaled Bondi power for different values of /3 and C,. 

luminosity for different values of C,, which describes the distribu- 
tion of the magnetic field in the space above the disc. It is found 
that the slopes of the model calculations with different values of 
are almost constant, which is in good consistent with the observed 
correlation (see Fig.[T](. The slope changes very little with differ- 
ent values of /3 (see Fig. This means that the slope of the ob- 
served correlation between Lboi/iEdd and ikin/iEdd is always 
consistent with our model calculations independent of the values 
of parameters adopted. Figure [3] shows that a relatively high C, is 
required for low-/3, which means that field lines diverse rapidly if 
the field strength is relative high, in order to explain the observed 
correlation between A and Lkin/iEdd- This provides useful clues 
to constructing detailed accretion disc/outflow models, which is be- 
yond the scope of this work. 

We find that the internal energy is not important compared 
with other terms in Eq. l |15t . The jet power is mainly related to 
the strength of the magnetic fields threading the disc and then to 
the disc properties. The bolometric power of the accretion flow is 
also related with the disc properties (e.g., density and temperature). 
Thus, the observed correlation between the kinetic power of the 
outflow and the bolometric power of the accretion flow can be nat- 
urally reproduced by our model calculations. The solid line in Fig. 
|4] corresponds to the best fitted model calculation on the observed 
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correlation between Eddington-scaled kinetic power and bolomet- 
ric luminosity, which implies that over 60% of the accreted gas at 
the outer radius escapes from the accretion disc in a wind before 
the gas falls into the black holes in these low luminosity AGNs. We 
find that best fitted model parameters C, and /? are somewhat degen- 
erated, i.e., they cannot be uniquely determined from the compar- 
ison with the observation, which may be caused by the simplified 
outflow model adopted in this work. The observed correlation can 
be used to constrain the model if the physics of the origin and the 
configuration of the magnetic fields is included, which is beyond 
the scope of this work. 

The results of our ADAF model with outflows can be qual- 
itatively explained in the frame of the self-similar ADAF model 
( lNaravan&"^l 19951) . The radiative efficiency o f an ADAF varies 
with mass accretion rate. iNaravan & yH ( Il995h 's results showed 
that the bolometric luminosity of an ADAF Ladaf oc m", where 
s ~ 2 for an ADAF due to the effects of energy advection in the 
accretion flow. Based on their self-similar solution for an ADAF, 
B'^ oc pgas cx: m and then Lj oc m, which leads to ij oc ^adaf ^ 

rO.5 

-iJADAF- 

Instead of ADAFs discussed above, one may wonder whether 
the standard thin accretion disc model can explain this correla- 
tion between A and Lkin/I/Edd, though these sources are accreting 
at very low rates. The dependence of jet power on dimensionless 
mass accretion rate for standard accretion discs with different scal- 
ing laws for magnetic field strength o n the disc p roperties were 
explored in some previous works (e.g., Moderski & Sikora 199^; 

ICa^ 



Ghosh & Abramowicj[T997l ; iLivio. Ogilvie. & Pringlelll999l : 



2002h. Most of them suggested that Lkin /I/Edd oc A (e.g., see 



Livio. Ogilvie. & Pringlell 19991 ; ICaoll2002l for the details), and no 



theoretical model calculations can reproduce the observed correla- 
tion: Lkin/j^Edd oc A° based on thin accretion disc m odels (e.g., 
iModerski & Sikordll99d iGhosh & Abramowicz|[l997l) . 

Our calculations also show that the kinetic power increases 
with Bondi power when we set the mass accretion rate at the outer 
radius A/o = i\/bondi (see Fig.O. This is roughly consistent with 
the correlation b e tween Pbondi/^Edd and -LkinZ-f'Edd found in 
iMerloni & Hein j ( |2007|) (see Fig. 1 in their paper). The mechan- 
ical efficiency e (Lkin = eAfbondiC ) varies from 10"^ to 10"^ 
in the AGN feedback simulations according to the complicated en- 
vironments and outflow s/jets mechanism, etc jOaspari et a"i]|2009l : 
iHopkins & slvislboiol) . In our calculations, the value of e is about 
10"^ ~ 10"^ which may be useful for future AGN feedback sim- 
ulations. 
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